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Using two extremely different models of glass formers in two and three dimensions we demonstrate 
how to encode the subtle changes in the geometric rearrangement of particles during the scenario 
of the glass transition. We construct a statistical mechanical description that is able to explain and 
predict the geometric rearrangement, the temperature dependent thermodynamic functions and the 
Q-relaxation time within the measured temperature range and beyond. The theory is based on 
an up-scaling to proper variables (quasi-species) which is validated using a simple criterion. Once 
constructed, the theory provides an accurate predictive tool for quantities like the specific heat or the 
entropy at temperatures that cannot be reached by measurements. In addition, the theory identifies 
a rapidly increasing typical length scale ^ as the temperature decreases. This growing spatial length 
scale determines the a-relaxation time as Tc, ~ exp{fj,^/T) where is a typical chemical potential 
per unit length. 
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I. INTRODUCTION 



For typical glass formers with soft potentials the "glass 
transition" is accompanied by a dramatic increase of re- 
laxation times: the viscosity of the super-cooled liquid 
shoots up by fifteen orders of magnitude within a rela- 
tively short temperature range |l|. When the liquid is 
examined on microscopic time scales, the molecules seem 
jammed, trapped into moving only within their cages. Of 
course, on the much longer viscous time scale these same 
systems do reach equilibrium by restructuring themselves 
through an ergodic exchange of particles between cages. 
As long as the system is ergodic it lends itself, at least in 
principle [2I], to a statistical mechanical treatment. Nev- 
ertheless, it is extremely hard to reach a useful statisti- 
cal mechanical theory as long as one considers the glass 
former on the level of its constituent particles because 
it is difficult to evaluate partition-function integrals in 
continuous coordinates. It is therefore tempting to find 
a reasonable up-scaling (coarse-graining) method that 
would define a discrete statistical-mechanics with par- 
tition sums rather than integrals, with the sum running 
on a finite number of quasi-species having well charac- 
terized degeneracies and enthalpies. Indeed, in a number 
of examples in 2-dimensions (2D) 0, 0, H, H, 0| and in 
one example in 3D [^] it was shown recently that such 
a discrete statistical-mechanics is possible and quite ad- 
vantageous d, in providing a successful description 
of the statistics and the dynamics of systems undergo- 
ing the glass transition. In this paper we examine how 
far the predictions of such a statistical mechanical theory 
can be pushed. As the theory is developed on the level of 
up-scaled quasi-species, one can ask whether the energy, 
entropy and other thermodynamic quantities appearing 
in the quasi-species language are identical numerically to 
the corresponding quantities of the system itself when 
computed on the level of constituting particles. One of 
the aims of this paper is to answer this question in the 
affirmative. This means that the statistical mechanical 



theory can be used to quantitatively predict physical ob- 
servables of interest also beyond the range of measure- 
ments, hopefully providing a satisfactory understanding 
of the phenomenology of glass formation. 

To demonstrate the generality of the approach we ex- 
amine two extremely different models of glass formers 
to which we apply the same procedure of: i) identify- 
ing an appropriate up-scaling ii) validating that the cho- 
sen up-scaling yields a consistent statistical mechanics 
and iii) demonstrating how the resulting theory provides 
a good understanding of the subtle geometric organiza- 
tional changes as a function of temperature, of the ther- 
modynamic functions and of the relaxation time. We 
will stress the point that understanding the way that the 
concentrations of the available up-scaled quasi-species de- 
pend on the temperature is tantamount to understanding 
the scenario of the glass transition. To underline the gen- 
erality of the approach we deal below with two models, 
one model is studied in 3D in an NPT ensemble and the 
other in 2D in an NVT ensemble. A short announcement 
of the results pertaining to the first model was published 
in [8] and here we provide considerably richer details. 
The results of the second model are novel to the present 
paper. The models were chosen, as explained below, to 
represent extremely different microscopic characteristics; 
nevertheless once up-scaled, the statistical mechanics ap- 
pears very similar encouraging us to propose that the 
approach is rather general. 

The structure of the paper is as follows: in Sect. |TT] 
we introduce the two models studied below, discuss the 
dynamics of their correlation functions as a function of 
temperature, and measure the relaxation time Tq that 
will be later connected to the statistical mechanical the- 
ory. In Sect, mil we present the upscaling for these two 
models. This upscaling defines quasi-species in terms 
of carefully chosen groups of particles rather than indi- 
vidual particles. We then show how the scenario of the 
glass formation can be characterized by the temperature 
dependence of the concentrations of the various avail- 
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able quasi-species. Since this upscaling is not unique, we 
present a validation of the choice of quasi-species. The 
vahdation is a demonstration that the choice of quasi- 
species provides a self-consistent statistical mechanics. 
In other words, each quasi-species is characterized by its 
own enthalpy and its own degeneracy. The statistical 
mechanics is then shown to be able to accurately recap- 
ture the scenario of the glass formation in temperatures 
where data is available, and to continue to predict it for 
temperatures that are outside the range of numerical sim- 
ulations. In Sect. IIVI we relate the upscalcd picture to 
the dynamics of the system. We show how a natural 
static length-scale is appearing, and how this length scale 
characterizes the relaxation time Tq. Finally, in Sec. |V] 
we demonstrate that the statistical mechanics as defined 
on the upscalcd quasi-species is able to provide the cor- 
rect thermodynamic function for the original system. In 
other words, we can compute the energy, entropy, spe- 
cific heat etc. directly from the upscalcd picture, in the 
range of temperatures accessible to simulations as well as 
for temperatures that are inaccessible to simulations. In 
Sect. IVII we discuss the degree of generality of the present 
approach and stress what are the remaining riddles for 
future research. 



II. TWO MODELS OF GLASS FORMATION 

A. Three-dimensional Binary Model in NPT 
ensemble 

Some results concerning the upscaling of the first 
model that is studied below were announced already in 
and here we provide further details. The model is 
a version of a much studied model [H, [H, [13, [lB| , here 
of a 50:50 mixture of N point-particles in 3-dimensions 
(N — 4096 in our case), interacting via a binary poten- 
tial. We refer to half the particles as 'small' and half as 
'large'; they interact via the potential U{rij): 



ao 







< 



> 



rdij) 
rc{ij) 



(1) 

Here, e is the energy scale and ct^ = l.Oa, 1.2(t or 1.4cr 
for small-small, small-large or large-large interactions, re- 
spectively. For the sake of numerical speed the potential 
is cut-off smoothly at a distance, denoted as Tc, which 

is calculated by solving dU / drij\r^-=r^ = which trans- 
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lates to Tc — {a/P)"-!^'' Uij. The parameter ag is cho- 



sen to guarantee the condition U{rc) = 0. Below we 
use a = 8 and /? = 6, resulting in Tc = ^^8/6 ct^ and 
ao = 0.10546875, see Fig. [TJ Note that the potential is 
purely repulsive, in a premeditated distinction from the 
second model discussed in Subsect. IIIBI 

The model dynamics were studied as a function of the 
temperature keeping the pressure fixed at p = 10 (NPT 
ensemble). The model was simulated by employing the 




FIG. 1: Color online: The pair-wise potential of the binary 
model 
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FIG. 2: Color online: Time dependence of the correlation 
functions ([2]) for a range of temperatures (decreasing from 
left to right) as shown in the figure. The inset shows the 
relaxation time Tq in a log-lin plot vs l/T, compared to an 
Arrhenius temperature dependence. 



Verlet integration scheme with the Berendsen thermostat 
and barostat [lH . The units of mass are the mass of the 
particles m, energy is measured in units of e and the 
fundamental length scale is a. The slowing down in the 
super-cooled regime is exemplified by measu ring the self- 
part of the intermediate scattering function [14| summed 
over the large particles only. 
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In Fig. [2] we show these correlation functions for k — 
5.1(7"^ and for a range of temperatures as indicated in 
the figure. We see the usual rapid slowing down that 
can be measured by introducing the typical time scale Tq, 
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FIG. 3: The pair-wise potential for the hump model. In the 
inset we show a snapshot of the position of the point-particles 
(the circles represent the radius a) . Note that the two typical 
scales rniin and r* appear as the typical distances between 
particles 



that is determined by the time where Fk {t — ; T) = 
Fk{0; T)/e= 1/e. The relaxation times are shown in the 
inset of Fig. [5] as a function of 1/T in a log-lin plot to 
stress the non-Arrhenius dependence at lower T. 



B. Two-dimensional 'hump' model in NVT 
ensemble 



The second model is constructed following [l6| such 
as to have very different microscopic properties from the 
binary model. The interaction potential is constructed 
as a piecewise function consisting of the repulsive part 
of a standard 12-6 Lennard- Jones potential connected at 
to a polynomial interaction P{x) = '^^QiX^. 
The flj's are tuned so that P{x) displays a peak at 
f — ''hump and also such that there is a smooth conti- 
nuity (up to second derivatives) with the Lennard-Jones 
interaction at C/(ro) = eft-o as well as with the cut-off in- 
teraction range [/(j'*) = 0. The interaction potential for 
the hump model reads: 



3 + ^0 , < ro 

ro < Tij < 
, rij > Ti, 

(3) 

with the parameters as shown in table 1. 

Note that the two typical distances that are defined 
by this potential, i.e. the distance at the minimum rmin 
and the cutoff scale r^,, appear explicitly in the amor- 
phous arrangement of the particles in the supercooled 
liquid, as shown in the inset in Fig. [31 The model has 
two crystalline ground states, one at high pressure with 
a hexagonal lattice and a lattice constant of the order of 
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TABLE I: Parameters used in the potential of the hump model. 



Tniin. At low prcssurc the ground state is a more open 
structure in which the distance r^, appears periodically. 
At intermediate pressures the system fails to crystallize 
and forms a glass upon cooling. 

The great difference between these two models will be- 
come even clearer after we define the up-scaled quasi- 
species below. On the face of it, the scenario for the 
glass transition appears similar as exemplified in Fig. [4] 
in which the intermediate scattering function 
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is shown for k — 6.16(t^^. The hump model was sim- 
ulated using the same algorithms as the binary model, 
but without the Berendsen barostat that was unneeded 
for the NVT ensemble. 

The relaxation time Tq for the hump model is shown 
in Fig. m in which we stress again the difference be- 
tween the region in which there is Arrhenius behavior 
and the region where the relaxation time grows faster 
than exp{E/T). 



III. CHOICE OF UP-SCALED VARIABLES AND 
VALIDATION 

Our aim is to provide a theory that captures, when 
the temperature goes down, the subtle changes in the 
structural organization of the particles; this rearrange- 
ment is in fact directly responsible for the glass transi- 
tion. The first step in our approach is up-scaling, (or 
coarse-graining) from particles to quasi-species that can 
be characterized by their enthalpy and degeneracy. Up- 
scaling is an art, since it can be done in various ways and 
there is no unique algorithm to select a-priori a 'best' 
up-scaling. Therefore we need to validate the choice of 
up-scaled variables using a criterion that was introduced 
in Q. The up-scaling is typically different for different 
models, and we describe it separately for the two models 
at hand. 
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Small Particles 
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FIG. 4: Color online: Time dependence of the correlation 
function [4] in the hump model for a range of temperatures as 
shown in the figure. Inset: the temperature dependence of 
the relaxation time Tc,{T). 
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FIG. 5: Color online: Temperature dependence of the con- 
centrations of the quasi-species of the binary model. Symbols 
are simulation data and the lines are a guide to the eye. 



A. Up-scaling of the binary model 



The potential in the case of the binary model is purely 
repulsive and it has well defined range of interaction 
I'cihj) for any given pair of particles. A natural up- 
scaling for this potential is provided by the particles and 
their nearest neighbors, where 'neighbors' are defined as 
all the particles j around a chosen central particle i that 
are within the range of interaction rc{i,j). We refer to 
the type of a central particle (small or large) in combi- 
nation with the amount of its nearest neighbors to define 
the quasi-species. In the interesting range of tempera- 
tures we find 8 quasi-species with one 'small' central par- 
ticle and 3, 4... 10 neighbors, and 9 quasi-species with 
one 'large' central particle with 6, 7 ... 14 neighbors, all 
in all 17 quasi-species. Other combinations have negligi- 
ble concentration (< 0.5%) throughout the temperature 
range. We denote these quasi-species as Cs{n) and Cein) 
with s and £ denoting the small or large central particle, 
while n denotes the number of neighbors. We measured 
the mole-fractions {Cs{n)){T) and {Ci{n)){T) and the 
results are shown in Fig. [5l We see that some concentra- 
tions decrease upon decreasing the temperature, others 
increase, and yet some first increase and further decrease. 
We submit to the reader that the subtle changes in con- 
figurational arrangement as seen in this plot encode the 
scenario of the glass transition in a way that we will at- 
tempt to decode below. 



B. Up-scaling of the hump model 



The potential of the hump model has a distinct min- 
imum which suggests that quasi-species could consist of 
each particle and all its neighbors which are within the 
minimum. This type of up-scaling may have a limited 
usefulness when the temperature is sufficiently high to 
allow easy escape over the maximum. The glass tran- 
sition occurs however in a temperature range where T 
is lower than the maximum, and we expect the chosen 
up-scaling to be sensible throughout the interesting tem- 
perature range. In this range we find quasi-species with 
2, 3, 4, 5 and 6 neighbors within the minimum. In Fig. [5] 
we show the temperature dependence of the quasi-species 
of this model which are denoted as Ct with i = 2, 3, . . . , 6 
being the number of neighbors within the minimum. The 
reader should be sensitive to the distinction between the 
two models. The quasi-species in the binary model can 
change upon every cage vibration; anytime the distance 
between two particles exceeds or goes below Tc the def- 
inition of the quasi-speices changes. On the other hand 
in the hump model the quasi-species are much more sta- 
ble, since an exchange of a particle calls for an escape or 
penetration by going over the hump. We have chosen the 
models to have these clear differences in order to test the 
generality of our approach. 
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FIG. 6: Color online: Temperature dependence of the con- 
centrations of the quasi-species of the hump model. Symbols 
are simulation data and the lines are a guide to the eye. 



C. Validation of the up-scaling 

Since there is no unique algorithm to choose the up- 
scaling, we need to have a criterion that validates or re- 
jects the choice of quasi-species. In ref. Q such a crite- 
rion was introduced as explained next. 



1. Validation of the up-scaling in the binary model 

To decide whether the up-scaling provides a useful 
statistical mechanics for the binary model we now ask 
whether there exist free energies JFs(n;T) and Ti(n]T) 
such that 



{Cs{n))(T) = 
{Ci{n))iT) = 



-rs(n;T)/T 



2y' 



-y^s{n;T)/T 



(5) 



The free energies are found by inverting Eqs. ([5|) in terms 
of the measured concentrations. In doing so one can al- 
ways choose one of the concentration to have by defi- 
nition zero free energy. This is because the constraint 
X]„C's('^) + SnC'f("') — 1 makes the system of equa- 
tions overdetermined. Then the free energies of all the 
other quasi-species are actually the difference from this 
chosen reference particle. We then plot these quantities 
as a function of the temperature, as demonstrated for the 
present case in Fig. [7l We say that our upscaling is vali- 
dated if J- sin; T) and Tiln; T) can be well approximated 
as linear in the temperature. Then we can interpret 
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FIG. 7: Color online: The approximate linear dependence of 
the free energies of the chosen quasi-species on the temper- 
ature. From the slope we read the degeneracy and from the 
intercept the enthalpies (up to normalization), cf. Eq. |6l 
Note that when the free energies are large we do not have 
data: the concentrations become exponentially small and in 
a finite simulation box they disappear completely. 



where now the degeneracies gs{n) and gi{n) (read 
from the slopes in Fig. [7]) and enthalpies Ti,s{n) and 
He{n) (read from the intercepts) are temperature- 
independent. This validates the choice of up-scaling. 
In other words, the approximate linearity of the inverted 
free energies in the temperature means that we can write 
the concentrations as 



{Cs{n)){T) 
{Ci{n)){T) 



g,(n)e-'^°(")/^ 
ff,(n)e-^^(")/^ 



(7) 



Then we can use these forms also as a prediction for 
temperatures where the simulation time is too short to 
observe the relaxation to equilibrium. The resulting de- 
generacies gs{n) and gi{n) can be easily modeled by a 
Gaussian distribution around the most probable number 
rimp of nearest neighbors for small and large particles 
respectively: 



9s{n) 
9e{n) 



3-[("-<p)V2T,^] 
= -[("-»Lp)V2'T?] 



n^p = 4.65, erf = 1.55 



7.50, cr| = 2.0 .(8) 



(6) 



The analytic fit to the measured degeneracies is shown 
in Fig. [HI upper panel. The same figure shows in the 
middle panel the enthalpies of the various quasi-species. 
One could model the enthalpies as a linear function in 
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FIG. 8: Color online: Upper panel: The degeneracies gs{n) 
and ge{n) read from the slopes of Fig. [7] (in circles) and the 
degeneracies according to the gaussian model Eq. Mid- 
dle panel: the measured enthalpies. Lower panel: compari- 
son of the measured concentrations of quasi-species to those 
calculated from Eqs. ([7} using the model degeneracies and 
measured enthalpies. Here symbols are data and lines are 
theoretical predictions. 



n. These results are easily interpreted; we have high 
enthalpies when there are large free volumes (few neigh- 
bors). The lowest enthalpies are found when there are 
many neighbors and there is no much costly free volume. 
In other words, at the present density and range of tem- 
peratures the pV term dominates the energy in the en- 
thalpy. Using the theoretical degeneracies and the mea- 
sured enthalpies we compute the concentrations of all our 
quasi-species and compare them with the measurement in 
the lowest panel of Fig. [H The agreement that we have, 
especially considering the number of quasi-species and 
the simplicity of the theory, is very satisfactory. Notice 
that the competition between degeneracy and enthalpy 
explains the rather intricate temperature-dependence of 
the concentrations of the quasi-species, sometimes declin- 
ing when the temperature drops, sometime rising, and 
sometime having non-monotonic behavior. 



2. Validation of the up-scaling in the hump model 

Validating the up-scaling in the hump model follows 
the same ideas as in the binary model. We choose the 
quasi-species n = 2 as the reference concentration with 
F2 = and then invert the data for {Ci){T) to find the 




FIG. 9: Color online: The approximate linear dependence of 
the free energies of the chosen quasi-species on the tempera- 
ture in the hump model. From the slope we read the degener- 
acy and from the intercept the energies(up to normalization), 
cf. Eq. [lOl 
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TABLE II: Temperature independent enthalpies and degen- 
eracies of the sub-species of the hump model. 



free energies via 



(9) 



The observed linearity of J'i{T) in the temperature means 
that we can write 



(10) 



where Tii is the enthalpy and gi the degeneracy of the 
ith quasi-specie. The test of linearity from which we can 
read the energies and degeneracies is shown in Fig. (5] 
From the intercepts of the lines in Fig. [2] we read the 
enthalpies, and from the slopes we read the degeneracies. 
The resulting numerical values are shown in Table |TT1 

Since the hump model is studied in the NVT ensem- 
ble the reader may ask why we get enthalpies rather 
than energies. The reason is that the upscaling using 
the index i takes into account only the particles that 
reside inside the hump. The energy of a quasi-species 
with i particles within the hump is proportional to i. 
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and thus for low temperatures the energetically preferred 
subspecies would be the one with no neighbors within the 
hump. Clearly this arrangement does not satisfy the con- 
stant volume constraint which needs to be taken into ac- 
count. Indeed, what we call enthalpies cannot be straight 
energies, otherwise the energies should have been pro- 
portional to i, making the smallest rather than the 
largest. 



D. Summary of the statistical mechanics 

We can summarize the findings up to now by saying 
that at least in terms of capturing the scenario of the 
changes in the spatial organization of the particles our 
statistical mechanics appears very suitable. Our quasi- 
species appear to have very well defined enthalpies and 
degeneracies and we can therefore capture the temper- 
ature dependence of the concentrations of the quasi- 
particles with high accuracy. Our thesis is that this sce- 
nario is actually the scenario of the glass transition, and 
hidden in it is also the reason for the slowing down, as 
demonstrated in the next section. Intuitively the picture 
should be clear already now. As the temperature re- 
duces the quasi-species with high free energy tend to dis- 
appear, leaving us with quasi-species of low free energy. 
Then the dynamics become more constrained, since it is 
more costly to change objects of low free energy (creat- 
ing on the way objects of higher free energy) than at high 
temperature when there are many objects with high free 
energy that can readily change to objects of lower free 
energy. What we need now is to make this observation 
more quantitative. 



IV. RELATION TO DYNAMICS 

In this section we connect the structural theory to the 
dynamical slowing-down. To this aim we note that in 
both models there are a number of quasi-species whose 
concentration goes down exponentially (or maybe faster) 
when the temperature decreases, and that the relaxation 
time shoots up at the same temperature range. We re- 
fer to these quasi-species as the 'liquid' ones; These are 
the quasi-species that are prevalent when the tempera- 
ture is high but are getting rare when the temperature 
goes down. We discuss now the connection between this 
phenomenon and the slowing down. 




FIG. 10: Color online: The temperature dependence of 
Ciiq(r) is shown as the upper continuous line. The contribu- 
tions of the various 'liquid' sub-species are shown with sym- 
bols which are identified in the inset. 



scale. 



e(T) ^ [pCu^in 



-1/3 



(11) 



where p is the number density. This length scale has 
the physical interpretation of the average distance be- 
tween the 'liquid' quasi-species. It was argued before 
0) IE 111 that this length scale can be also interpreted as 
the linear size of relaxation events which include 0(^(T)) 
quasi-species. We can therefore estimate the growing free 
energy per relaxation event as AG = P-S,{T) where p is 
the typical chemical potential per involved quasi-species. 
This estimate, in turn, determines the relaxation time as 



(12) 



The quality of this prediction can be gleaned from 
Fig. [TTl where we can see that the fit is excellent, with 
fi « 0.3. The intercept in Fig. [TT]is of the order of unity; 
this is very reassuring, since this is what we expect when 
T — > c». 



B. The hump model 



A. The binary model 

In this example the 'liquid' concentrations are those 
consisting of small particles with 3-7 neighbors, and large 
particles with 6-11 neighbors, see Fig. [10] [ll|. We sum 
up these concentrations and denote the sum as (Ciiq)(r). 
The dependence of (Ciiq)(T) on the temperature is shown 
in Fig. [ini This concentration is used to define a typical 



In this model the quasi-species whose concentration 
goes to zero rapidly in the relevant temperature range 
are those with i = 2, 5 and 6. In Fig. [T^lwe see their tem- 
perature dependence and also their sum which is again 
denoted as Ciiq(r). Since in this model we are in two 
dimensions, the typical scale £,{T) is related to Ciiq(T) 
according to 

^ [pCuq(T)]-l/2 ; (13) 
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FIG. 11: Color online: The relaxation time Ta{T) in terms 
of the typical scale ^{T) in the binary model. We show the 
excellent fit to Eq. (I12|l with /i — 0.37. Note that the intercept 
at T cxD is of the order of unity as it must be. 
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FIG. 12: Color online: The quasi-species contributing to 
Ciiq(r) in the hump model, and the temperature dependence 
of the their sum, which is Ciiq(r). 

With this obvious change we expect Eq. to be valid 
also here, and indeed in Fig. [13] we see that this expec- 
tation is wonderfully fulfilled. 

C. Pertinent Remarks 

A few points should be stressed. As we expect (cf. Ref 
0), in systems with point particles and soft potential, 
there is no reason to fit the relaxation time to a Vogel- 
Fulcher form [l[ which predicts a singularity at finite tem- 
perature. In our approach we predict that ^ —>■ oo only 
when T ^ 0, and there is nothing singular on the way. 
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FIG. 13: Color online: The relaxation time Ta{T) in terms 
of the typical scale (,{T) in the hump model. We show the 
excellent fit to Eq. ((13} with fi = 0.34. Note that the intercept 
at T ^ oo is of the order of unity as it must be. 



only slower and slower relaxation. At some point the 
simulation time will be too short for the system to re- 
lax to equilibrium, but we can use Eq. (|12p to predict 
what should be the simulation time to allow the system 
to reach equilibrium. It is important to bear in mind 
that finite systems of point particles with soft poten- 
tial are different from finite granular media or systems 
of hard spheres which can truly jam and lose ergodic- 
ity. Point particles with soft potential remain ergodic 0] , 
and therefore should be amenable in their super-cooled 
regime to statistical mechanics. We argue that in or- 
der to construct simple, workable statistical mechanics 
one needs to up-scale the system and find a collection of 
quasi-species with well defined enthalpies and degenera- 
cies. In this paper we showed that the simple criterion 
to validate the choice of up-scaling which was suggested 
in P can be applied to very different models once we se- 
lect the proper up-scaling. Once the structural theory is 
under control, a natural length scale appears and can be 
used to determine the relaxation time, also for tempera- 
tures that cannot be simulated due to the fast growth of 
the relaxation time. The fact that the present approach 
works equally well in two and three dimensions, in NPT 
and NVT ensembles and in very different models provides 
good reason to believe that it has a substantial degree of 
generality. 



V. THE THERMODYNAMICS OF GLASS 
FORMING SYSTEMS 

As a last issue we raise the following question: is the 
statistical mechanical theory as defined on the up-scaled 
quasi-species determining also the thermodynamics of 
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FIG. 14: Color online: The ratio of the LHS of Eq. (fT4)l . 
denoted as AHth + C (here, C = 2.8) over the RHS of Eq. 
(|14[) . denoted Affcxp. We conclude that the enthalpy of the 
system is excellently reproduced by the up-scaled variables. 



FIG. 15: Color online: Comparison of the LHS of Eq. ([TSjl in 
red continuous line, to the RHS of the same equation, in sym- 
bols. We conclude that the energy of the system is excellently 
reproduced by the up-scaled variables down to temperatures 
of about T « 0.05 For lower temperatures the extrapolation 
obviously fails. 



the original system on the particle level. The answer 
to this question is not obvious a-priori, but turns out to 
be in the affirmative. As done throughout this paper, we 
demonstrate this point separately for the two models at 
hand. 



A. Thermodynamics of the binary model 

The binary model is studied in the NPT ensemble and 
therefore the first question to ask is whether the average 
enthalpy of the quasi-species represents the correct en- 
thalpy of the system. In other words, we have extracted 
the enthalpies Hgiji, T) and Hi{n, T). Is it then true that 
the enthalpy of the system, computed as the average po- 
tential energy summed over all the pairs of particles and 
summed with pV should equal the enthalpy of the quasi- 
species: 

10 14 

Y,Cs{n)ns{n,T ) +J2Ci{n)nfXn,T) 

n— 3 n— 6 

In answering this question we remember that our en- 
thalpies are defined up to an arbitrary constant since 
we chose one of the quasi-species as a reference with zero 
free energy. We can therefore add or subtract an arbi- 
trary constant from the LHS or the RHS of Eq. [TH In 
Fig. [H we show the ratio of the LHS of Eq. ^ to the 
RHS, with the conclusion that the enthalpy of the sys- 
tem is excellently reproduced by the up-scaled variables. 
Needless to say, we could therefore compute the specific 
heat at constant pressure, Cp, from either set of data, 
and from that we could get the entropy of the system, 



reassuring us that the up-scaling method, once validated 
as above, provides also the correct thermodynamics for 
the glass forming system. In the next subsection we show 
actual computations of such thermodynamic quantities. 



B. Thermodynamics of the hump model 

The hump model was studied in the NVT ensemble, 
and therefore it is interesting to test whether one obtains 
the correct system's energy. The potential energies of the 
quasi-species which were characterized by the number of 
neighbors that are bound within the minimum arc to 
an excellent approximation linear in that number, Ei — 
0.83i as measured with respect to the same zero point 
as that of the potential We therefore need to check 
whether 

ooo 6 N 

— J2ain+T = — J2{u{n,)) . (15) 

We have added T for the excitation around the ground 
state of the quasi-species (two degrees of freedom in 2D). 
This test is shown in Fig. [Ij] where the predictions of 
the theory were extrapolated all the way to zero temper- 
ature. Obviously the data loses its predictability some- 
where around T — 0.05 where the specific heat becomes 
negative. This artifact is seen even better in the plot of 
the specific heat which is offered in Fig. [161 Nevertheless 
the plot of Ct, vs. T indicates the expected typical spe- 
cific heat peak around T = 0.1 which could not be seen 
in the simulations due to the inherent limitations of com- 
puter time. Once the specific heat takes a dive, at some 
point we can no longer believe the extrapolation, and 
more accurate data are necessary to be able to predict 
the thermodynamic properties near zero temperatures. 
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FIG. 16: Color online: Symbols: the derivative with respect 
to T of the RHS in Eq. (|15|l . Continuous line: derivative with 
respect to T of the LHS in Eq. (|15|l . The theory predicts the 
expected specific heat peak and fails when the specific heat 
becomes negative. 

VI. SUMMARY AND THE ROAD AHEAD 

In summary, we have shown how the scenario of the 
glass transition can be neatly characterized by the tem- 
perature dependence of the concentrations of the various 
quasi-species that are obtained after up-scaling. Those 
quasi-species that are readily changed because they are 
rich in free energy are of course those that deplete first. 
We are left with quasi-species that are low in free energy, 
and these would naturally be hard to change since their 



change necessarily increases the free energy. The pre- 
sumed ergodicity allows us to introduce statistical me- 
chanics which can beautifully encode the scenario, us- 
ing a small number of quasi-species together with their 
enthalpies and degeneracies. The statistical mechanics 
applies also at temperatures that are beyond the range 
available to molecular simulations, and can be used to 
predict what is happening there. We have shown that 
when the upscaling and the statistical mechanics work 
well, the thermodynamics of the system can be under- 
stood very well using the up-scaled picture. 

There are still issues to understand. One of our main 
concerns is how to achieve a first-principles prediction of 
the parameter /i in Eq. ((T2)) . We have, so far, been un- 
able to bridge this gap, leaving the connection between 
the dynamics and the upscaling theory still dependent 
on a phenomenological fit. A satisfactory prediction of /i 
would, in our opinion, close the loop and put to rest the 
riddle of the slowing down in the glass- forming systems. 
Also, the upscaling approach had been so far limited to 
computer model of glass formation, we propose that it 
would be enlightening to try and apply it to real exper- 
imental systems. We hope that these and other issues 
would be clarified in future research. 
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